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SUMMARY 

High-speed helicopter rotor impulsive noise prediction is an important problem of aeroacoustics. The 
deterministic quadrupoles have been shown to contribute significantly to high-speed impulsive (HSI) noise 
of rotors, particularly when the phenomenon of delocalization occurs. At high rotor-tip speeds, some of 
the quadrupole sources lie outside the sonic circle and move at supersonic speed. Brentner has given a 
formulation suitable for efficient prediction of quadrupole noise inside the sonic circle. In this paper, we 
give a simple formulation based on the acoustic analogy that is valid for both subsonic and supersonic 
quadrupole noise prediction. Like the formulation of Brentner, the model is exact for an observer in the 
far field and in the rotor plane and is approximate elsewhere. We give the full analytic derivation of 
this formulation in the paper. We present the method of implementation on a computer for supersonic 
quadrupoles using marching cubes for constructing the influence surface (E surface) of an observer space- 
time variable (x, t). We then present several examples of noise prediction for both subsonic and supersonic 
quadrupoles. It is shown that in the case of transonic flow over rotor blades, the inclusion of the supersonic 
quadrupoles improves the prediction of the acoustic pressure signature. We show the equivalence of the new 
formulation to that of Brentner for subsonic quadrupoles. It is shown that the regions of high quadrupole 
source strength are primarily produced by the shock surface and the flow over the leading edge of the rotor. 

The primary role of the supersonic quadrupoles is to increase the width of a strong acoustic signal. 

1. INTRODUCTION 

High-speed flight of helicopters has been an illusive goal because of the rapid increase in rotor noise and 
vibration. Noise and vibration increase dramatically because of the increasing disparity between the fluid 
velocity over the advancing and the retreating rotor blades. When the supersonic how region on the advancing 
rotor blade extends off the blade into the far held (a phenomenon known as delocalization [1]) Thigh-speed 
impulsive (HSI) noise becomes dominant over all the other rotor noise sources. For this reasonTthe efficient 
prediction of the HSI noise of helicopter rotors is currently an important problem of aeroacoustics. The cause 
of this noise has been identihed since the late 1970’s [112] as the deterministic quadrupoles in the vicinity 
of the rotorTandrin the case of delocalizationT beyond the sonic circle and the blade tip. Many schemes 
have been proposed by researchers based on the acoustic analogy [3-9] and the Kirchhoff method [1011 1]. 
At presentTmost of these schemes are limited to subsonic quadrupole source motionrhoweverTvery recently 
several methods have been proposed for the prediction of noise from supersonic sonic sources [12-15]. In 
the case of subsonic quadrupole noise predictionTa method based on a formulation by BrentnerTformulation 
QlATexists [7] which is exact for an in-planeTfar-held observer but is approximate elsewhere. This method 
has been implemented in the code WOP WOP + [617] and has been shown to be highly efficient and robust. 
There is a need for an efficient and robust method of prediction of the supersonic quadrupole noise. We 
present a new formulation — based on the same model used by Brentner — that is valid for both subsonic and 
supersonic quadrupole noise prediction. 

The derivation of a supersonic quadrupole formulationFwhich we call formulation Q2ris the main result of 
this paper. In Section 2 we start with an exact solution of the wave equation for quadrupole sources of the 
Ffowcs Williams-Hawkings (FW-H) equation given by Farassat and Brentner [16]. In this solutionFvolume 
integrals involving the quadrupole sources are only differentiated with respect to observer time. We write 
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these volume integrals in terms of a surface integral over the collapsing sphere and a source time integral. For 
an observer in the rotor plane and in the far fieldrthe collapsing sphere is approximated as a right circular 
cylinder normal to the rotor disc and the quadrupole source strength is integrated along lines normal to 
the rotor disc and treated as sources on the rotor disc. We then hypothesize that the quadrupole noise 
everywhere can be predicted using these surface sources. This hypothesis has been validated by Brentner 
and Holland [6] . The idea of approximating the volume (quadrupole) sources with equivalent surface sources 
was originally proposed and numerically implemented by Yu et al. [1] for the far-held solution of the FW-H 
equation. Later this idea was also implemented by Schultz and Splettstoesser [3]TBrentner and Holland [6 ]T 
and Brentner [7]. Our main contribution has been to use this idea in obtaining closed-form solutions of the 
same equationrformulations Q1A and Q2Fwhich seem to result in more efficient prediction of HSI noise. The 
new formulation (Q2) presented here is very simple and unlike formulation Q1A is valid for both subsonic and 
supersonic quadrupoles. Hence we can directly compare quadrupole noise predictions from both formulations 
in the subsonic case. 

In Section 3 we discuss how formulation Q2 is implemented in a new testbed code called WOPWOP2+. 
The quadrupole sources beyond the sonic circle can have multiple emission times and the usual solution of 
the wave equation for subsonic surface sources (e.g.rformulations lriArand Q1A) will have a singularity 
known as the Doppler singularity. To avoid this Doppler singularityT it is necessary to use a E-surface 
formulation [17118] . The method of construction of the E surface used in WOPWOP2+ is known as marching 
cubes [19] — an algorithm originally developed for computer graphics [20]. 

Some examples of HSI noise prediction for a hovering rotor is presented in Section 4. First we present a 
comparison of the results of noise prediction for subsonic quadrupoles based on formulations Q1A and Q2. It 
is shown that the results agree well with each other. A study of the surface source strength (i.e.Tthe integral 
of quadrupole source strength along the line normal to rotor disc) shows that the primary contributions to 
the quadrupole noise come from the shock surfaces on and beyond the blade and the flow over the leading 
edge of the blade. It is also shown that the inclusion of the quadrupole sources beyond the sonic circle 
improves the prediction of the width of the main pulse and the shape of the acoustic pressure signature and 
agreement with experimental data. FinallyFwe demonstrate the robustness of the formulation by performing 
predictions for out-of-plane and near-field observers. Concluding remarks follow in Section 5. 

2. FORMULATION AND SOLUTION OF THE PROBLEM 

We begin with the solution of the following wave equation for quadrupole noise radiation from the Ffowcs 
Williams-Hawkings equation: 


= d^ [TijH{f)] (1) 

where the bar over the partial derivative operator indicates generalized different iationFT^ is the Lighthill 
stress tensorr H (/) is the Heaviside functionr and f — 0 describes the blade surface (/ > 0 outside the 
blade). The solution for this equation was given by Farassat and Brentner [16] as follows: 


4t rp' Q (x,t) 


+ 


+ 



-oo /> 0 


— ^ dQ dr 
r 


Off 3 T rr - T u 

at J J 

-oo /> 0 

c J J 3Trr ~ Tii dtl dr 

-oo /> 0 


( 2 ) 


where the quantity T rr is the double contraction TijfifjT and fi are the components of the unit vector in the 
radiation direction. In additionrdQ is an element of the collapsing sphere surface g — 0. We now assume 
that the observer is in the far field and on the rotor plane. The part of the collapsing sphere intersecting 
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Figure 1. The actual and approximate collapsing sphere surfaces in the vicinity of the rotor blade. 


the source region near the rotor blade can be approximated by a right circular cylinder normal to the rotor 
plane. This is shown in figure 1. 

Let us assume that the rotor is nominally in the ^/i^-plane and y% isT thereforeTperpendicular to this 
plane (i.e.Tthe rotor tip-path plane). We integrate the inner integrals of equation (2) with respect to y% over 
the approximate collapsing sphere surface. Let us define 


Qij(yuV 2 ,T) = / Tij dy 3 
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and use the relation 


dT = c drdY (4) 

where dT is the length of an element of the curve defined by the intersection of the collapsing sphere with 
the rotor disc. Equation (2) can now be written as 

Id 2 f 1 
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where the subscript ret indicates the integrand is evaluated at the retarded time. The integrals in equation (5) 
are all over the entire 2 /i?/ 2 -plane. This fact allows us to bring the observer-time derivatives inside the integrals 
without worrying about the limits of integration. 

The next step is the most crucial in the derivation of the final result. We note that d/dt = d/dt | x 
meaning x in the frame fixed to the undisturbed medium is kept fixed in this differentiation. Also note that 
the tensor Q with components Qij is defined in the x frame of referenceThence all velocity terms in Qij(x, t) 
are specified with respect to the frame fixed to the undisturbed medium. HoweverTif we use a change of 
coordinates (x,£) — >• (? 7 ,r)rwhere the 77 -frame is aligned with the rotor blade (i.e.Tby coordinate rotation 
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and translation) Twe have 
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where 77 is the position vector in the rotating frame and r is the source time. Here V = dr] /dr is the 
velocity of the point 77 specified in the frame fixed to the undisturbed medium. We note that V is in the 
rotor plane. It is important to recognize that when we refer to Qij\ v we really mean that the components 
of the tensor Qij represented in coordinates that are instantaneously aligned with the rotating frame. Thus 
equation (6) provides the time derivative of Qij | x in the stationary frame in terms of Qij\ v which is specified 
in the coordinates of the moving frame. Using this operator notation in equation (5)Twe get 
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Notice that the operator L r operates on Qij only because Ti and fj do not depend upon t or r. To write 
this equation in final form we express V(?7,r) as follows: 

V = V /.' + a; x 77 (8) 

where Vf is the forward velocity of the rotorTo; is the angular velocity of the rotorTand 77 is the position 
vector of the source Qij in the rotor plane with the origin at the rotor center. We assume that both V p and 
< jJ are time independent and we note that u = |u;|. From equation (8)Twe have 


ThusTwe can express L 2 T as follows: 
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All the gradients are with respect to rj. We have 
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where d/dr) is the directional derivative in the 77 (radial direction) and 77 = |t 7|. We also have 
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and V = (Vi, V 2 ). ThereforeT !/ 2 can be written as 
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When this expression is used in equation (7)Twe get a singularity free expression for supersonic quadrupole 
noise prediction. We note that equation (7) has second space and time derivatives of Qij as well as first 
space derivatives in the /yi^-plane (the rotor plane). These quantities are available in the CFD postprocessor 
that is used to compute Qij for acoustic calculations. We will refer to equation (7) as formulation Q 2 . As 
it standsT formulation Q 2 is valid for subsonic and supersonic quadrupole noise prediction for helicopter 
rotors in hover or forward flight. Note that this equation is very simple and has no singularities. We have 
assumed that the shocks on the blades are smeared over one or more grid cells and that Qij has continuous 
second derivatives over the rotor plane (although the magnitude of the second derivative can be very high at 
the foot of the shocks). These assumptions are generally satisfied in CFD calculations for helicopter rotor 
aerodynamics. 

2.1. ANALYSIS OF MAIN RESULT 

We will now do an order of magnitude study of the far-held term of our main resuhTequation (7). We can 
draw very useful conclusions from such a study as will be shown below. We know that the peak of directivity 
of quadrupole noise is in the rotor plane with the observer ahead of the helicopter. Let us put the observer 
in such a location in the far held. Then in the frame hxed to the undisturbed medium (where aq-axis is the 
hight direction and the £ 3 -axis is normal to the rotor plane) T the components of the unit radiation vector 
can be approximated as r = (1,0,0). The major contribution to the far-held quadrupole noise comes from 
Q 11 which we will look at closely below. 

The numerator of the integrand of the far-held term is 
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We can now estimate the order of magnitude of each term in equation (15) as follows. Let the advancing tip 
speed be denoted by Vat • ThenTwe see that 
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In these equationsTthe derivative djdr\\ is the directional derivative in the chordwise direction. Note that 
for a hovering rotor dQn/ drYd 2 Qn/ dr 2 YdJvd Vf are all zeroTtherefore the only remaining component is 
(V • V') 2 Qii. Since uYin generalTis small for helicopter rotorsTwe can see that the dominant term in forward 
hight is also most likely the term (V • V / ) 2 Qn. The right side of equation (16c) can be further estimated as 


,2 d 2 Q n V 2 t Q n 

^ dnl ~ (A ^) 2 


(17) 


where A/q is the chordwise scale over which signihcant change in Qn occurs. Signihcant changes of Qn 
occur near both the leading edge stagnation point and the base of the shock. Hence we suspect that the 
dominant sources of quadrupole noise will be located at these locations (A/q ~ LE radius and A/q ~ width 
of the projection of the shock surface in the rotor plane). These conclusions are verified by the computation 
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Figure 2. Contours of L^Qn in the vicinity of a hovering UH-1H rotoiT M // = 0.925. 

of L 2 T Qn for a hovering UH-1H model rotor blade at tip Mach number 0.925 shown in figure 2. It is 
apparent in the figure that the primary source of HSI noise is the shock wave (as proposed by Farassat and 
his colleagues — see references 16121-23) and the flow over the leading edge of the blade. The significance of 
the quadrupole source in the leading edge region has not been widely recognized in previous workThoweverT 
it has been found here to be of comparable amplitude to the quadrupole source strength in the shock region. 


3. NUMERICAL IMPLEMENTATION 

A new codeT called WOPWOP2+ris used to demonstrate the utility of formulation Q2. WOPWOP2+ 
differs significantly from WOP WOP + [617] in that it uses a E-surface formulation to compute thickness and 
loading noiseFas well as the quadrupole noise. The construction of the E surface and subsequent integration 
over the E surface is performed using the method of marching cubes integration developed by Brentner [19]. 
The numerical calculation of quadrupole noise has been divided into two stages: a preprocessing stage in 
which the integration of the Lighthill stress tensor in the normal directionT indicated in equation (3)Tis 
carried outTand an evaluation stage in which the quadrupole contribution to the acoustic pressure specified 
in equation (7) is determined. Both the preprocessor and the acoustic calculation are described briefly in this 
section. More information on the preprocessorTwhich is the same preprocessor that is used by WOPWOP+T 
can be found in reference 6. 

3.1. PREPROCESSOR 

Although the evaluation of Qij can be performed independently of observer position and retarded timeTthe 
preprocessor must read in the CFD solutionrinterpolate the solution at the necessary quadrature locationsr 
and then perform the numerical quadrature in the direction normal to the rotor disc. The preprocessor needs 
knowledge of both the CFD grid topology and the solution format. In the implementation used for this workr 
the interpolation of the CFD data is two dimensional and is done one radial station at a time. For a given 
radial stationrdata are interpolated to quadrature points needed for composite Gauss-Legendre integrationr 
on lines normal to the rotor plane. The lines are uniformly distributed in the chordwise direction. A two- 
dimensional linear least-squares interpolation is used to interpolate the densityTmomentumrand energy at 
each quadrature point. The Lighthill stress tensor Tij is evaluated with the interpolated data. The value of 
Qij on the rotor disc is determined at each chordwise location before moving to the next radial station. The 
results are stored for the acoustic calculation stage. The CFD data used for the noise prediction [24] was 
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Figure 3. Typical quadrupole grid used for WOPWOP2+ calculations for hovering UH-1H rotor. Note that 
only every 4th grid line is shown in the chordwise direction. 

given on a grid with rather coarse radial resolution beyond the blade tip. The quadrupole source strength 
radial resolution off the blade tip was increased five fold by utilizing 6th order polynominal interpolation 
along radial grid lines for both the location and strength of the quadrupole data. A typical grid generated 
by the quadrupole preprocessor is shown in figure 3. 

3.2. WOPWOP2+ 

The primary function of the WOPWOP2+ code is to perform the integration indicated in equation (7) 
numerically. Although the integration is over the entire rotor disc planeTin practice the source strength 
is zero over a large part of the plane; henceTthe quadrupole integration is only performed near the rotor 
blade (see figure 3). The main difficulty in the numerical evaluation of equation (7) is the construction of 
the E surface for the portion of the plane in which the source strength is nonzero. The E surface is the 
collection of points in space-time that emit signals that reach the observer at one particular observer time. 
The integration is complex because the pointwise mapping between the physical source plane and the E 
surface is not known explicitly. Special care must be taken in the construction because in practice the E 
surface may be composed of several disjoint pieces when the source motion is supersonic — exactly the case 
we are interested in. Figure 4 shows a typical evolution of the E surface. Notice that in in the first few time 
steps each position on the E surface is single valued (i.e.Ta point only occurs at one azimuthal position)r 
but later the leading edge and eventually the entire outer portion of the surface becomes triple valued (i.e.T 
point near the tip of the grid may occur at three azimuthal positions). 

The marching cubes method begins constructing the E surface by choosing the source time and computing 
the corresponding observer time and integrand value at each grid point. If the observer times are computed 
and stored for each desired source timeTthe discrete computational data become a three-dimensional array; 
two computational indices parameterize the surface spatially and a third index accounts for the source time. 
In this three-dimensional computational spacer isosurfaces of observer time t areT by definitionT distinct 
realizations of the E surface. The extension of the marching-cubes algorithm for surface integration [19] 
determines how the surface intersects a logical cube in the three-dimensional computational gridTcomputes 
the contribution to the integral from that portion of the surfaceTand then moves (or marches) to the next 
cube. The topology of the surface within a single cube can be determined uniquely by examining the 
function value (observer time in this case) at each of the cube vertices and comparing this value to the 
desired surface value. A table lookup is then used to determine the exact topology of the surface in the 
current cube. The surface is formed by a set of triangular panels that have vertices on the edges of the 
cubes. The value of the surface integral over each triangle is approximated as the average integrand value of 
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Figure 4. Evolution of the E-surface geometry and integrand strength distribution TiTjL 2 T Qijjr a function 
of observer time. The observer is located 3.09R from the rotor hub and the rotor tip velocity is Mh = 0.95. 
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the triangle vertices multiplied by the triangle area. Linear interpolation is used to determine the integrand 
values at the triangle vertices based on the previously computed value at the cube vertices. The marching 
cubes algorithm is generic — the only difference in computing thicknessrioadingTor quadrupole noise is the 
value of the integrand computed at each grid vertex. (For more detail on the marching cubes algorithmFsee 
references 19 and 20.) 

The version of the marching-cubes algorithm used for the calculations in this paper was adaptive — that 
is the integral value for a single cube was compared to the value obtained using the 8 subcubes obtained by 
bisecting each of the sides of the original cube. If the difference between the integral values was greater than 
a specified toleranceTeach of the subcubes were also subdivided in a recursive manner. Six levels of recursion 
were allowed in the computations for this paper. When a cube was dividedTthe quadrupole source strength 
at the new vertices was computed using linear interpolationTbut the observer time and source location were 
computed without approximation. Finally the integrand values were computed anew at each of the new 
vertex locations. 

For simplicityTthe current WOPWOP2+ code only implements formulation Q2 for a hovering rotor. For 
a hovering rotorFL r Qij and L\Qij can be written in the form 

LrQij = ~(u xt}) ■ V n Qij = (18a) 

Of]i 


and 


L 2 T Qij = (w X n) • V„((w X T}) • v„Qy) = ( 18b ) 

where djdr\\ is the directional derivative in the azimuthal direction and r\ — \r]\. Equation (18) is implemented 
numerically in WOPWOP2+ with a second-order accurate central-difference operators. 

4. NUMERICAL RESULTS 

In this section we first make a comparison of the acoustic pressure signatures of a hovering model rotor 
from WOPWOP+ and WOPWOP2+ which use formulations Q1A and (irrespectively. ThenFwe present 
the comparison of the predicted and measured acoustic pressure signatures for the same hovering model 
rotor at four different tip Mach numbers. At the three highest tip Mach numbersr the phenomenon of 
delocalization occurs and WOPWOP2+ must be used for HSI noise prediction if the contribution of the 
supersonic quadrupoles is to be included in the prediction. The remaining analysis examines the role of the 
supersonic quadrupoles and the robustness of the formulation. 

A model-scale rotor test conducted by Boxwell et al. in 1978 [25] and repeated later by Purcell in 1988 [26] 
is used for comparison. The measured data was for nonlifting hovering rotor generating HSI noise. The rotor 
was a l/7th scale UH-1H main rotor with straight untwisted blades and NACA 0012 airfoil section. The 
rotor radius R was 1.045 m with a chord of 7.62 cm. The measured data reported here are all from a 
microphone in the rotor plane and at 3.09 R from the rotor tip. The Euler solutions utilized as input in this 
numerical work were provided by Baeder and are described in references 5 and 24. The Euler solutions are 
also used for direct comparison with the acoustic prediction when experimental data is unavailable. 

4.1. COMPARISON OF FORMULATIONS 

Figure 5 shows a comparison of the predicted acoustic pressure signatures from WOPWOP+ (formulation 
Q1A) and WOPWOP2+ (formulation Q2) at tip Mach number 0.925. The thickness and loadingFquadrupoler 
and total acoustic pressure time histories predicted by each of the codes are also shown in this figure. 
Although the supersonic quadrupoles are important in prediction of the acoustic pressure signature because 
delocalization occurs at this operating conditionrWOPWOP+ can only handle subsonic quadrupole sources; 
thereforeTwe have only used the subsonic quadrupole sources in both predictions for this comparison. We 
have utilized the marching cubes approach to construct the E surface (influence surface) of the rotor blade 
and the quadrupole source surface in WOPWOP2+. A good agreement in this comparison proves two 
points. Firstrit will tell us that the construction of the E surface is correct in WOPWOP2+. Secondrthe 
two formulations Q1A and Q2 are equivalent. Both these points are evident in figure 5. This figure also 
shows that the individual components due to thickness and loading and the quadrupole sources as well as the 
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Figure 5. Comparison of WOPWOP+ prediction (•) with subsonic WOPWOP2+ prediction ( ) for a 

UH-1H hovering rotorPM# = 0.925. Quadrupole grid extends 0.075 R beyond rotor tip. (a) thickness and 
loading components; (b) quadrupole component; (c) total acoustic pressure. 
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total acoustic pressure signatures from the two codes agree well. ThusFwe have established some confidence 
in using WOPWOP2+ for prediction of HSI noise. 

4.2. COMPARISON WITH MEASURED DATA 

We now present HSI noise calculations for tip Mach numbers 0.88r0.9r0.925Fand 0.95 in figure 6. The 
quadrupole grid extends 1.86 R beyond the blade tip for all the WOPWOP2+ calculations shown in figure 6. 
For comparisonrwe have also shown the WOP WOP + signature which includes quadrupole sources up to the 
sonic circle. It is seen that the agreement of the WOPWOP2+ signature with the measured data is excellent 
and better than that of WOPWOP+ for each case. We haverthusrdemonstrated the ability to predict the 
noise from supersonic quadrupoles in the case of delocalized shocks and the resulting improvements in the 
overall shape and level of the acoustic pressure. 

For the more intense cases ( Mh > 0.90)rthe agreement of the WOPWOP2+ prediction with the mea- 
sured acoustic pressure signature is not fully satisfactory because the WOPWOP2+ prediction overpredicts 
the negative peak pressure. This is apparent in figure 6 for the times between the WOPWOP+ and the 
WOPWOP2+ shock locations. A closer examination of the E surface and the quadrupole integrand strength 
distribution for this time rangeTshown in figure 7Freveals that this time range is precisely when the region of 
the strongest quadrupole integrand strength is is bifurcating (i.e. transitioning from subsonic to supersonic). 
Notice in particular that the area of the highest and lowest values (lightest and darkestTrespectively) change 
dramatically during this time range. Hence relatively small errors in the peak values of L^Qij can result in 
amplified error. To reduce the sensitivity of the acoustic predictions to this type of errorTwe utilized a 7 
pointTmoving-window-average smoothing in the chordwise direction of the quadrupole source strength. This 
smoothing had only a small effect on the acoustic signal away from the shockrbut was useful to reduce the 
error previously noted. Figure 8 is a representative comparison of the smoothed and unsmoothed acoustic 
predictions. (Note smoothing was used in all of the computations for this paper except the Mh = 0.88 case 
where the shock strength did not seem to warrant smoothing.) 

It is well known that the quadrupole accounts for the nonlinear propagation effects caused by the finite 
particle velocity and the variation of sound speed in the physical problem. The calculations shown in figure 
6 seems to indicate that the primary role of the supersonic quadrupoles is to increase the pulse width of 
intense propagating waves. The width of the main pulse of the signatures predicted by WOPWOP2+ agrees 
very well with the measured signature. Even for the Mh = 0.88 caseTwhich is not delocalizedrthe supersonic 
quadrupoles improve the agreement by increasing the width and amplitude of the acoustic signal. For the 
delocalized casesrthe supersonic quadrupoles also decrease the slopeFand thus improve the agreementrof the 
triangular shape proceeding the rapid increase in acoustic pressure. Figure 9 shows the effect of the extent 
of the quadrupole grid for the particularly intense Mh = 0.95 case. Three separate computations were 
made with the quadrupole source grid extending beyond the rotor tip 0.05itT0.33itTand 1.86iHrespectively. 
The quadrupole grid is shown in figure 9a with the three grid extents indicated. The first WOPWOP2+ 
computationr0.05i? beyond the rotor tipTis essentially identical to the WOPWOP+ calculation shown in 
figure 6d. Notice that the prediction for the largest grid extentn.86i? beyond the rotor tipTagrees very well 
with the data in both waveform width and shape. These computations suggest two things: even though the 
quadrupole correctly predicts the nonlinear propagation (i.e.Tthe widening and changing of the waveform 
shape) Tit is probably an inefficient tool for predicting the nonlinear propagation because an accurate CFD 
computation must proceed the quadrupole prediction; and it seems likely that the nonlinear propagation 
could be determined more appropriately by another method starting with the acoustic signal somewhat 
closer to the rotor. 

4.3. FORMULATION ROBUSTNESS 

In this section we wish to demonstrate the robustness of the method by performing predictions which 
violate some of the assumptions leading to formulation Q2. First we predict the noise for out-of-plane 
observer and compare the acoustic pressure with an Euler solution. SecondlyFwe predict the noise for a very 
near-field observer. 

The acoustic pressure was predicted using WOPWOP+ and WOPWOP2+ at two observer locations 
directly below the an in-plane observer at 3.09 R from the rotor hub. These observers are 10 and 20 deg below 
the rotor planeFrespectively. In figure 10 we show these prediction for the Mh = 0.9 case. No measured data 
are available at these observer locationsMiereforerwe have interpolated the Euler solutions [5] used as input. 
It is seen that the two sets of calculations agree fairly well with each otherTbut the WOPWOP2+ result 
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Figure 6. Comparison of WOPWOP+ ( ) and WOPWOP2+ ( ) predicted acoustic pressure with 

experimental data [26] (□) for hovering model UH-1H rotor. Quadrupole grid in WOPWOP+ prediction 
extended almost to sonic circle and in WOPWOP2+ predictions extended 1.86/? beyond the rotor tip. (a) 
M h = 0.88; (b) M h = 0.90; (c) M H = 0.925; (d) M H = 0.95. 
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Figure 7. Evolution of quadrupole integrand strength distribution fiTjL 2 T Qijjr on the emission surface as a 
function of observer time for times near the peak negative acoustic pressure. Dark and light shading indicates 
negative and positive values of rirjL^Qij/r respectively. Note the extremely rapid change in the contours 
area as the source region bifurcates. (a) t = 1.050 msec; (b) t = 1.117 msec; (c) t = 1.183 msec; (d) t = 1.250 
msec. 
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Figure 8. Comparison of WOPWOP2+ prediction with unsmoothed Qij input data (•) with the prediction 

made using Qij smoothing ( ). A 7-point moving- window average smoothing was applied in the chordwise 

direction for a UH-1H hovering rotorPM# = 0.925. 


characterizes the waveform width and asymmetry slightly better for the 10 deg down observer. Furthermorer 
as expectedrthe peak negative value of the main pulse of the acoustic pressure reduces with increasing 
observer angle. When the observers is 20 deg below the rotor plane the agreement is good for both methodsr 
but the quadrupole contribution is relatively small. The chordwise smoothing of the quadrupole source 
strength in the WOPWOP2+ computation is believed to be responsible for the lower amplitude in the 
10 deg down prediction; since the quadrupole contribution is small in the 20 deg down case the effect of 
smoothing is not readily noticeable. 

In a second comparisonFshown in figure HTwe have predicted the acoustic pressure with WOPWOP2+ 
at an in-plane observer 1.094 R from the rotor hub. This observer is just inside the sonic circle and is less 
than two chordlengths from the rotor tip at the closest distance. The quadrupole grid extends almost to 
the observer. WOPWOP+ was unable to determine the retarded time satisfactorily for this severe test case 
because the root finding algorithm exceeded the number of iterations allowed. The WOPWOP2+ prediction 
slightly underpredicts the Euler solution; neverthelessTthe comparison is extremely good at this very near- 
field location. In particularTthe WOPWOP2+ prediction agrees perfectly with the Euler data both before 
the negative peak and for the curved part of the signal at the top of the shock-like structure after the 
negative peak. The far-held and near-held quadrupole terms (terms with 1/r dependence and terms with 
1/r 2 and 1/r 3 dependenceFrespectively) from equation (7) are also shown in hgure 11. Clearly the near-held 
quadrupole terms — usually neglected by other researchers — contribute signihcantly to the correct prediction 
of waveform shape at this close distance. Figure 11 demonstrates the importance of keeping all of the terms 
so that the acoustic prediction can be compared directly with CFD [16]. Both hgures 10 and 11 demonstrate 
the robustness of formulation Q2. 


5. CONCLUDING REMARKS 

We have presented a new quadrupole noise prediction method based on a new analytic resultT called 
formulation Q2Tvalid for both subsonic and supersonic quadrupole sources. The new formulation is very 
simple and without any singularity. The procedure for implementation of the result is discussed in the 
paper. This new code is called WOPWOP2+. We have demonstrated that formulation Q2 is equivalent to 
formulation Q1A of Brentner used in WOPWOP+ for subsonic quadrupole sources. By order of magnitude 
study of the formulation Q2 far-held integrandTwe have shown that the shock surfaces and the stagnation 
how at the leading edge of the blade are regions of high source intensity. We have shown that for rotors 
operating at high tip Mach numbers — before and after delocalization — the new formulation predicts acoustic 
pressure signatures which agree well with the experimental data in both the shape and the level of the main 
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Figure 9. Effect of quadrupole grid extent shown for hovering UH-1H rotorYMn = 0.95. (a) quadrupole 
grid (every 10th grid line in chordwise direction shown); (b) predicted acoustic pressure for three quadrupole 

grid extents: T0.05R beyond rotor tip (subsonic); T0.33R beyond rotor tip; T1.S6R beyond 

rotor tip; □ experimental data [26]. 
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Figure 10. Comparison of WOPWOP+ ( )rWOPWOP2+ ( )rand Euler [5124] (•) predicted acoustic 

pressures for hovering model UH-1H rotorPM# = 0.90. (a) 10 deg below rotor plane; (b) 20 deg below rotor 
plane. 
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Figure 11. Near-field noise prediction at in-plane observer located 1.094 R from rotor hub for hovering model 

UH-1H rotorr Mh = 0.90 with quadrupole grid 0.0935 R beyond rotor tip. T far-held quadrupole 

component; r near-field quadrupole components; rWOPWOP2+ total acoustic pressure; •rEuler 

[5124] total acoustic pressure. 


pulse of the signature. We have also shown that the supersonic quadrupoles widen and modify the shape 
of the waveform. A new and robust option is now available for prediction of HSI noise of helicopter rotors 
based on formulation Q2. 
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